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Abstract 

We undertake an investigation of particle acceleration in the context of non-linear electrodynamics. We 
deduce the maximum energy that an electron can gain in a non-linear density wave in a magnetised plasma, 
and we show that an electron can ‘surf’ a sufficiently intense Born-Infeld electromagnetic plane wave and 
be strongly accelerated by the wave. The first result is valid for a large class of physically reasonable 
modifications of the linear Maxwell equations, whilst the second result exploits the special mathematical 
structure of Born-Infeld theory. 


1 Introduction 

The implications of theories that couple the electromagnetic field to itself have been an enduring source of 
interest to particle theorists for decades, and recent developments in ultra-high intensity lasers have led to a 
surge of interest in relativistic non-linear electrodynamics by the wider community. It is expected that facilities 
such as ELI Qj will permit investigation of laser-matter interactions at intensities where non-linear perturbative 
effects mediated by virtual electron-positron pairs will be evident [2] , and numerous studies of the implications 
of the Euler-Heisenberg Lagrangian have been undertaken in this context. In the longer term, it is anticipated 
that vacuum pair production, a fundamentally non-perturbative QED process, will become accessible in the 
laboratory (for a recent review of the attendant theoretical issues see Ref. ED- In addition to QED processes, 
attention has also been paid in recent years to tests of axion electrodynamics using ultra-high intensity lasers 0]. 

Considerable progress has been made during recent years in the exploitation of large-amplitude plasma 
waves for particle acceleration. Such schemes are particularly attractive in the laboratory because the electric 
fields in a plasma wave can be several orders of magnitude greater than those sustainable in standard radio¬ 
frequency accelerator cavities. The most prevalent schemes realised thus far employ the strong fields in the wake 
behind an intense laser pulse propagating through the plasma [5] , although electron-driven wakefields have also 
been exploited for electron acceleration [Bj. Furthermore, recent developments have focussed on proton-driven 
wakefields as a paradigm for efficiently accelerating leptons to TeV energies [7]. In the astrophysical context, 
plasma waves were recently invoked to explain the emission of energetic electrons from within the interiors of 
pulsars [5]; such electrons are necessary for the formation of the electron-positron plasma populating a pulsar’s 
magnetosphere. The magnetic fields found in neutron stars are typically ~ 10 8 T, whilst those in magnetars 
may be two orders of magnitude higher, and non-linear effects due to QED are expected to be significant in 
such environments [Uj. One may speculate that non-Standard Model couplings also play an important role. 

The Born-Infeld Lagrangian is probably the most famous theory of non-linear electrodynamics whose mo¬ 
tivation lies outside the Standard Model. It first appeared in the 1930s as a classical model of the electron 
with finite self-energy [TUI, resurfaced during the mid-1980s as an effective action in string theory m and 
its notoriety was finally cemented during subsequent years in the context of D-branes [12] - The Born-Infeld 
Lagrangian has a privileged mathematical status within the family of Lagrangians that depend only on the 
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electromagnetic field tensor, the dual of the electromagnetic field tensor and the spacetime metric tensor (but 
not their derivatives). It is the only such regular generalization of the vacuum Maxwell Lagrangian whose field 
equations exhibit an absence of birefringence and shocks mm- 

Such considerations led to a recent study nn of the properties of non-linear density waves in a Born-Infeld 
plasma and their implications for particle acceleration. An approximation to the maximum energy gain of an 
electron trapped and accelerated by the wave was obtained, but the result was independent of the Born-Infeld 
coupling parameter. Section [2] below revisits this topic from a general perspective with an investigation of 
particle acceleration in a non-linear density wave propagating along the ambient field lines of a magnetised 
plasma with arbitrary electromagnetic self-couplings. We show that the non-appearance of the Born-Infeld 
coupling parameter in the estimate of the maximum energy gain of the trapped electron in Ref. m has nothing 
to do with the exceptional properties of the Born-Infeld Lagrangian. Indeed, we demonstrate that an exact 
calculation of the maximum energy gain yields the same result for any physically reasonable Lagrangian that 
is an algebraic expression of the two fundamental invariants of the electromagnetic field, and this result is 
independent of the strength of the ambient magnetic field. 

Unlike many other non-linear theories, the source-free Born-Infeld field equations possess numerous non¬ 
trivial exact solutions which are indispensible for exploring non-perturbative aspects of the theory. In particular, 
an exact solution describing an electromagnetic pulse immersed in a uniform magnetic field is known m and 
we suggest that this may have implications for vacuum laser acceleration in future facilities, such as ELI [l|. 
Section [3] shows that an electron interacting with a non-linear Born-Infeld electromagnetic plane wave can be 
uniformly accelerated to arbitrarily high energies. This novel result is non-perturbative and has no analogue in 
linear Maxwell electromagnetics. 

Heaviside-Lorentz units are used throughout the following with c = 1 and, to avoid an unnecessary plethora 
of indices, intrinsic geometrical notation is used extensively. Further details of the notation and conventions 
used here may be found in Ref. [ 1 8j . 


2 Particle acceleration in a strongly magnetised plasma 

Much attention has been devoted to uncovering the behaviour of particles trapped in an electron density wave 
driven by an intense laser pulse, or particle bunch, propagating through a plasma. Although fully 3-dimensional 
configurations (the ‘bubble regime’) now pervade such studies [HIED]; the original laser-wakefield accelerator 
concept was formulated using 1-dimensional considerations |21] and approaches that employ non-linear plane 
waves remain useful for providing estimates. Furthermore, one can argue that plane waves are appropriate for 
modelling particle acceleration in neutron star crusts. The strong magnetic field within a neutron star polarizes 
the iron outer crust of the neutron star and leads to a highly anisotropic conductivity [22]. Electron density 
waves excited within the magnetic flux tubes are essentially free to propagate along the magnetic field, but 
their motion transverse to the held is greatly restricted; hence, to a first approximation, it is sufficient to only 
consider motion along the magnetic held lines 0. 

The purpose of this section is to show that the maximum gain in energy of an electron in a non-linear plane 
wave in a magnetised plasma is invariant within a large class of theories of electromagnetism. The plasma 
is modelled as a superposition of two charged pressureless perfect fluids, where one fluid describes mobile 
electrons and the other describes the charge carriers of a neutralizing background medium. It is assumed that 
the spacetime curvature is negligible and the worldlines of the charge carriers of the neutralizing background 
are timelike geodesics. Hence, we adopt the Minkowski metric 

g = —dt (g) dt + dx ® dx + dy ® dy + dz ® dz (1) 

and choose the 4-velocity Vo of the neutralizing background to be Vo = dt- 

The latter approximations are reasonable if the motion of the background is negligible over the timescales of 
interest. This is certainly the case if the background consists of ions whose charge-to-mass ratio is approximately 
three orders of magnitude lower than that of the mobile electrons (as in a wakefield accelerator), or if the 
background is the polarised outer crust of a neutron star and non-inertial terms (due to the rotation of the 
neutron star) in the metric components can be ignored. 
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The worldlines of the mobile electrons are trajectories of the unit normalised future-pointing timelike 4- 
vector field V. It is assumed that the motion of the electrons is parallel to a constant ambient magnetic field 
and their 4-acceleration X/yV satisfies 

VyV = —l v F (2) 

m 


with 


g{v,v) = - i (3) 

where qiyF is the Lorentz 4-force acting on the mobile electrons, iy is the interior product with respect to 
V, —q = e is the elementary charge, m is the electron rest mass, F is the electromagnetic 2-form and V is 
the Levi-Civita connection. The 1-form V is the metric dual of the vector field V i.e. the 1-form V satisfies 
V ( U ) = g(V, U ) for all vector fields U. Although the electromagnetic self-force (radiation reaction) on the 
electrons plays a significant role in strong fields [23], little is known about the self-force outside of the context 
of classical linear electromagnetism or perturbative QED. The emphasis of the following is on maintaining 
generality in the electromagnetic sector, and so the Lorentz force is adopted in the absence of a theory of 
radiation reaction in the general context. 

The background fluid is described by the electric 4-current enoVo where the constant no has the physical 
dimensions of a number density. In the case of a plasma, no/Z is the proper number density of the ions where 
Z is the ionisation multiplicity. The electromagnetic field equations may be written covariantly as 

dF = 0, d*G =—qn*V + qno*Vo (4) 

where n is the proper number density of the electron fluid and the Hodge map * is induced from the 4-form *1 
given as 

*1 = dt A dx A dy A dz. (5) 

A constitutive relation specifying the excitation 2-form G in terms of the 2-form F must be given to close the 
system of field equations ©, ([3]) , © . Equations ©, © , © may be derived [lB] from a Lagrangian containing 
a 0-form-valued function £em of the invariants X and Y defined as 


I = *(FA*L), Y = *(F A F), 
where 


*G = 2 


<9£em „ <9£em 

-ar* F+ sr F ) 


and the choice £em = X/2 yields classical linear Maxwell theory. 
For present purposes it is advantageous to replace © with 

dr K = jo A l k F 


( 6 ) 

(7) 

( 8 ) 


where jo = qn o * Vq and K is a Killing vector. The stress (-energy-momentum) 3-form tk is the sum of 
contributions from the electromagnetic field and the electron fluid: 

tk = LkF A *G + £em * K + mng(K , V) * V. (9) 

S ---' S -V-' 

EM field electron fluid 

Equation © expresses the local balance of energy, momentum or angular momentum when K is chosen appro¬ 
priately. For example, if AT is a unit timelike Killing vector then © connects the rate of change of the total 
energy density of the electromagnetic field and the electron fluid with the rate of work done per unit volume 
by the Lorentz force on the background fluid. Balance laws associated with linear or angular momentum are 
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revealed when K is chosen to be a generator of translations or rotations, respectively. It can be shown that the 
two systems of equations ©, ©, © and ©, ©, © are equivalent. 

The above model is applicable to a wave in the electron fluid density propagating along the lines of an 
ambient magnetic field whose curvature and evolution can be neglected. In this case, the electromagnetic field 
F is a superposition of an ambient homogenous static magnetic field (0,0, H) and the electric field (0,0, E{Q) 
driven by the electron fluid: 

F = -E(C) dt A dz — B dx A dy (10) 


where C, = z — vt is the phase of the density wave and the phase velocity v is a constant satisfying 0 < v < 1. 
Although no generality is lost by requiring v > 0, the choice |u| < 1 is a critical ingredient in the wakefield 
accelerator paradigm since particles must have the opportunity to be trapped in the wave and accelerated. 

It has long been known that non-linearities arising purely from the matter content lead to an upper bound 
on the amplitude of the electric field of a steady density wave whose wave 4-vector is spacelike [24}[25] . The 
maximum amplitude, known as the ‘wave-breaking limit’ by the plasma accelerator community, is an important 
parameter in the wakefield accelerator paradigm and is sensitive to the details of the plasma model muni. 
Although the magnetic field considered here is aligned with the velocity of the electron fluid and cannot directly 
influence their motion, it may influence their motion indirectly through electromagnetic self-coupling in the 
excitation 2 -form © and thereby affect the wave-breaking limit. 

It is convenient to analyse the field equations ©, ©, © using the pair {e 1 , e 2 }: 

e 1 = vdz — dt, e 2 = dz — vdt = d( (11) 

where the orthonormal coframe {'ye 1 ,'ye 2 ,dx,dy}, with 7 = l/\/l — v 2 , is adapted to observers moving at 
velocity v along z in the rest frame of the neutralizing background medium (i.e observers at rest in the ‘wave 
frame’). We seek a 4-velocity field V of the form 

V = K C) el +^(C)e 2 ( 12 ) 


where, using g(V, V) = —1, the component xjj is 

if) = — \//i 2 - y 2 


(13) 


with the sign of ip chosen to ensure that the velocity je 2 (V) of the plasma electrons in the wave frame is 
non-positive. Hence, the speed of the plasma electrons in the frame of the neutralizing medium is less than the 
phase speed v of the wave. 

The components of G in the basis {dt, dx, dy, dz} depend only on ( and it follows that d( A d*G = 0. Hence, 
© yields the electron proper number density n as a function of y, 


n 0 v"f 2 


(14) 


allowing n to be eliminated from the stress 3-form tk ■ 

Analysis of the balance law © proceeds by choosing the 1-form K to be each member of the basis 
{e 1 , e 2 , dx, dy} in turn. Both sides of © vanish when K £ {dx,dy} (the momentum transferred between 
the electromagnetic field and electron fluid is along z only) whereas 


tk — mnovy e 1 A dx A dy 


(15) 


for K = e 1 , with ~ indicating equality modulo closed forms and where n has been eliminated using (1141) . 
Insertion of (1151) into © leads to 


m 1 dy 
Q 7 2 ' 


(16) 
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Equations ([ 6 ]), (flOl) lead to X = E 2 — B 2 , Y = 2 EB, which with K = e 2 yields 


t k ^ 


2£2 ^em _ 2 E B d£RM 


dX 


dY 


+ C 


EM 


mriQVyJy 2 - y 2 


e 1 A dx A dy, 


(17) 


and 


dy 2 i 
jo A lkF = mno—e Ae AdxAdy 
di 


where ed has been used. 

Hence, equations (fT7|) . (flSl) give 


(18) 


dc[ 2E W 


+ 2 ELB 


<9£em 

9F 


£em - mn 0 {vyjy 2 - q 2 



(19) 


The periodic solutions to (ED for y have certain properties. Clearly y > 7 due to the square root in ED 
and the periodic solution with the largest amplitude satisfies y{Ci) = 7 where £ = is a zero of dy/dC, (and 
hence, using (fTTH) . E((i) = 0 ). Thus, the largest amplitude solution to (fT!Jl) satisfies 


2 E 2 ™ + 2 EB - £em - mno{v\Jy 2 - y 2 - y) = mnol - £em|b=o- (20) 

The maximum value /i ma x of y may be determined by evaluating ED at the turning point £ = £n of y 
immediately after the turning point ( = Ci- Using (fTfil) and dy/dC\^ u = 0 in (l20l) yields 

Mmax = Ai(Cll) = 7 3 (! + V 2 ). (21) 


Thus far, we have only considered the mobile electrons that form the wave (i.e. those described by the 
electron fluid). We now turn to the behaviour of electrons captured by the density wave. However, electrons 
trapped in the density wave do not, in general, follow the worldlines of the electron fluid (since, in contrast to 
the outcome of ED, they may propagate faster than the wave) and, in general, elucidating the impact of the 
captured electrons requires intensive numerical computation. However, if the population of captured electrons is 
sufficiently small then their back-reaction on the wave can be neglected. In this case, the motion of the captured 
electrons is dictated by the total electromagnetic field of the electron fluid, the neutralizing background and the 
background magnetic field only; hence, the captured electrons are modelled as test particles. 

The electric field is static in the wave frame and it follows that the quantity 

Ay = y(Cn) ~ MCi) = 27 V (22) 


is proportional to potential difference and therefore proportional to the increase in the energy of a test electron 
that moves from £ = (j to £ = (n. The coefficient of proportionality is straightforward to obtain using the 
following elegant argument. 

The test electron’s worldline C satisfies the Lorentz equation 


mX7 qC = qi(jF (23) 

where the electron’s 4-velocity C satisfies g(C,C) = —1. The unit timelike Killing vector K = 7 (d t + vd z ) 
satisfies V qK = 0 and it may be shown that d[g(C, K)\/dr = (V qC)(K) where r is the electron’s proper time. 
Hence d[—mg{C , K)]/dr = gi^i^F where — mg(C,K ) is the energy of the electron in the wave frame, and 
it follows that the change in energy of the electron over the interval [ri,rn] is q f^ n LqLkF dr. The previous 
integral may be written in the covariant and parameterisation-independent form q J c lkF and the change in 
energy A E K , 

AE K = -mg{C,K)\™, (24) 
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of the test electron in the wave frame follows immediately: 

AE K = q [ i k F = /-% dC = -Ag, = 21717 V. (25) 

Jc J 7 aC 7 
Ci 

Previous estimates of the maximum energy gain [SUE] were obtained using quantities, such as the maximum 
electric field of the density wave, that require £em(^, V) to be specified explicitly. In the above, we bypassed 
such details and obtained a lower bound on the electron’s maximum energy gain independent of £em(A, Y). 
Although the profile of the density wave depends on the details of £em(A, Y). the potential difference between 
two adjacent nodes of the maximum amplitude electric Held is independent of £em(A, Y). If the maximum 
amplitude density wave is weaker than its counterpart in classical Maxwell electromagnetism, then the period of 
the wave must be longer to ensure that (12511 holds. The above result is immutable and explains why the estimate 
of the maximum energy gain in the unmagnetised Born-Infeld plasma wave given in Ref. m is independent of 
the Born-Infeld parameter. 

Unlike the results in the wave frame, the energy gain of the test electron in the frame of the neutralizing 
background depends on the electron’s initial conditions. As with the electron fluid, we will assume that the test 
electron’s motion is along the magnetic field lines only. 

The electron’s 4-velocity C may be written as C = 7 M (AT + uL) where K = 7 (d t + vd z ), L = j(d z + vd t ), 
7 u = 1/vT — u 2 and (0,0, it) is the 3-velocity of the electron in the wave frame. Thus, the change in energy 
A Eg t of the electron in the frame of the neutralizing background is 

A E dt = —mg(C, d t )\ ™ = m 7 [ 711(1 + ujjv) - 71(1 + uiu)] (26) 

where ui, rtu and 71 , 711 are the values of u and 7 „, respectively, at the points I, II on C. Since —g(C 7 K)\^ = 
7 ii — 7 i follows trivially from the definition of K we find 

7n - 7t = 27 V 2 (27) 

using (ESD, and ( 1251 ) . ( 1271 ) can be used to express A Eg t in terms of 71 and v. The electron begins at a node of 
the electric field and we choose 77 > 0; it can then be shown that dAEg t /dyi < 0 and it follows that the largest 
value of AEg t . for fixed v, is lim 7I _>i AEg t = 4m(7 3 — 7). An electron that starts at a node of the electric 
field with speed u = 0 + in the wave frame, and reaches the adjacent node, gains energy AEg t = 4to( 7 3 — 7) in 
the frame of the neutralizing background medium. Thus, the maximum change in energy of an electron that 
begins at rest in the frame of the background is 4m(7 3 — 7) + 771(7 — 1 ) = m(47 3 — 37 — 1 ). Although the 
latter result was previously derived in the context of classical linear electromagnetism m , using a very different 
approach to that presented here, our novel observation is that the same result holds for more general theories 
of electromagnetism. 

Although the maximum energy gain does not explicitly depend on the details of the Lagrangian £em(A, Y) 7 
it may implicitly depend on £em(A, Y) since it depends on the phase speed v of the density wave. In particular, 
if the density wave is driven by an ultra-strong laser pulse then v is related to the structure and behaviour of 
the pulse which, in turn, are influenced by the details of £em(A1, Y) and the total electric 4-current. 


3 Particle acceleration in a Born-Infeld plane wave 

Several directions for generalizing the previous analysis are possible, and all involve an electromagnetic 2- 
form F whose structure is more general than m■ An option accessible to non-perturbative analysis is to 
focus on particle acceleration in regions where the total electric 4-current can be neglected. Indeed, the on¬ 
going evolution of ultra-high-power laser technology has greatly increased the feasibility of directly accelerating 
matter in free space (so-called ‘vacuum laser acceleration’) and the implications of an effective self-coupling of 
the electromagnetic field are of considerable interest in this context. 

Before turning to non-linear electromagnetism, it is worth noting that plane wave solutions to the linear 
vacuum Maxwell equations are not normally consider to be useful for particle acceleration [28]. In particular, 
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suppose that a pulsed plane wave collides with an electron. After a sufficient time interval, the pulse must 
completely overtake the electron (regardless of the electron’s initial conditions) because, in classical linear 
Maxwell theory, plane waves propagate rectilinearly at speed c = 1 without dispersing. However, it is common 
to choose an electric field whose profile is an antisymmetric function of the plane wave’s phase when exploring 
the behaviour of an electron driven by an intense few-cycle laser pulse EM, and therefore the integral of the 
electric held over the length of the pulse vanishes. It follows that there is no net change in the energy of the 
electron if radiation reaction is neglected. 

To be more precise, consider the 2-form F = £(z — t) (dz — dt ) A dx where the once-differentiable function 
£ : R —> R has compact support. The electric and magnetic fields of the pulse are (—£, 0,0) and (0, —£, 0) 
respectively. If the worldline C : r 4 (f(r), x(t), y(r), z(r)) of the electron satisfies the Lorentz equation 

toV qC = qi,jF with g(C, C ) = — 1, and C is the particle’s 4-velocity with r the electron’s proper time, then it 
is straightforward to show that mt(rp) = mt(ri) if J^ 1 £(<f>) dcj) = 0 where </> = z — t and i = dt/dr. It follows 
that the difference between the initial and final values of the electron’s relativistic energy mt vanishes. Vacuum 
laser acceleration relies on using a tightly-focussed ultra-high-power laser pulse to directly accelerate electrons 
in free space |31| . and the effects of the pulse cannot be adequately captured by modelling it as a plane wave; it 
is vital to account for the non-trivial pointwise dependence of the electromagnetic field on x,y. However, there 
is no reason to conclude that plane waves are ineffective for particle acceleration if the effects of the quantum 
vacuum (of the Standard Model or otherwise) are manifest. Indeed, as we will now show, a non-perturbative 
analysis in this context yields new results that are inaccessible using perturbation theory. 

The non-linear generalization of vacuum Maxwell electromagnetism introduced by Born and Infeld [lOj is 
a privileged theory. It is the only theory generated by a Lagrangian of the two electromagnetic invariants 
X : Y that agrees with vacuum Maxwell theory in the weak-field regime and whose solutions do not exhibit 
birefringence and do not develop shocks [SHU]- The Born-Infeld field equations also emerge from string/M 
theory cam and this has ignited modern interest nWTlIMTl in it as an effective theory of electromagnetism 
in strong fields. 

The source-free Born-Infeld field equations are 


dF = 0, d * Gbi = 0, 


(28) 


where the excitation 2-form Gbi = 2(<9y£bi F — dyC-Bi * F) is generated from the following 0-form: 

Gbi (X, V) = 4 (! - y/l - k 2 V - ^Y 2 /4) . (29) 

hi 

The Born-Infeld parameter k controls the strength of the self-coupling of the electromagnetic field, and vacuum 
Maxwell theory is recovered in the limit k —> 0. 

A substantial number of exact solutions to (12511 have been discovered mmm despite the non-linear 
structure of Born-Infeld electromagnetism. In particular, the exact solution m 


F = £(z — vt) (dz — vdt) A dx — B x dy A dz — B y dz A dx — B z dx A dy + \£(z — vt)dt A dz (30) 

to (1281) describes an electromagnetic plane wave propagating through an ambient uniform magnetic field 
(B x ,B y , B z ) where 


k 2 B z B x v _ ll + n 2 B 2 

1 + h 2 B 2 ’ V ~ V 1 + k 2 B 2 ’ 


(31) 


and B = yjB 2 + B 2 + B 2 . The wave propagates along the 2 -axis with phase velocity v and its profile is 
encoded by the function £ : R —> R. The electric field measured by an inertial observer with 4-velocity d/dt 
is (E x = —v£ 1 E y = 0 ,E Z = \£) and, in general, the electric field of a Born-Infeld electromagnetic wave has 
a non-zero longitudinal component in addition to the usual transverse component found in vacuum Maxwell 
theory. Remarkably, the smooth function £ is essentially unconstrained; the only requirement is that it satisfies 
the bound 


1 - k 2 X - k 4 Y 2 / 4 > 0 


(32) 
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arising from the argument of the square root in (12911 . 

Inspection of m shows that the phase velocity v of the wave is equal to or less than the speed of light 
c = 1 in vacuum Maxwell theory. Thus, it is possible for a particle to propagate faster than the phase speed 
v of the wave described by (1301) . Although a perturbative analysis in k would simply lead to corrections to 
predictions of linear Maxwell theory, we expect Born-Infeld electrodynamics to have novel implications resulting 
from non-perturbative considerations. Indeed, as we will now show, a test electron driven by the wave can 
have constant 4-acceleration and this result has no analogue in vacuum Maxwell electrodynamics. 

As in Section [2l the worldline C : (£(t), x(t), y(r), z(t)) of the test electron is presumed to satisfy the 

Lorentz equation 

mV^C = qtfjF (33) 

and the normalization condition 

g(C,C) = -1 (34) 

where C = id/dt + xd/dx + yd/dy + zd/dz with i = dt/dr, x = dx/dr , y = dy/dr. z = dz/dr. 

The properties of the Born-Infeld plane wave solution (1501) suggest that it may be fruitful to ask whether 
solutions to (1351) . (155)) exist that have constant phase, i.e. z — vi = 0. This condition may be written as CC, = 0 
where C = z — vt as before, and an electron satisfying CC, = 0 may be envisaged as ‘surfing’ the wave. However, 
although we are free to choose CC\ T =o = 0, there is no guarantee that a solution to (1551) . (1551) exists such that 
CC vanishes at r > 0. However, as we will now show, such solutions do exist if the properties of the profile £ 
are appropriate. 

In the following, it is convenient to introduce the type-(l, 1) tensor F that satisfies F(a,U) = a(n/F) for 
all choices of 1-form a and vector U. Furthermore, to avoid an unnecessary plethora of brackets or indices, 
it is useful to denote tensor contraction by juxtaposition; in particular, all = a(U ), aFU = F{a,U) and 
FQ = F(—, d a ) ® G{dx a , —) where G is a type-(l, 1) tensor and the Einstein summation convention is used with 
{dx a } = {dt, dx, dy, dz}, {d a } = {d/dt, d/dx, d/dy , d/dz} and a = 0,1, 2,3. 

The constraint CC = 0 may be written as 


dCC = 0 (35) 

and a consequence of (1551) is that its derivative S/^idCC) along C must vanish. Thus, a second constraint 

dCFC = 0 (36) 

is generated from the first constraint m using (El and V(jdC — 0. 

Since inspection of (151)1) reveals that is proportional to dCC, it follows that vanishes as a conse¬ 
quence of (1551) . Hence, differentiation of the second constraint (1361) along C yields the third constraint 

dCF 2 C = 0 (37) 

and, likewise, differentiation of (1571) along C leads to the fourth constraint 

dCF 3 C = 0 (38) 

and so-on, where F n = n p = iF. However, the Cayley-Hamilton theorem can be invoked to write F n , for n > 4, 
in terms of a linear superposition of F 3 , F 2 , F, and the identity tensor, so no new constraints are generated 
for n > 3. 

A condition on the electromagnetic field is obtained when (1551) . (1551) . (1571) . (1551) are regarded as a linear 
system for the components of C. Then, a non-trivial C exists if and only if A = 0 where the 4-form A is 


A = dC A dCF A dCF 2 A dCF 3 . 


(39) 


It can be shown that the 0-form *A is a cubic polynomial in £ with one simple root £ = 0 and one repeated 
root £ = £ cr it where the critical value f C rit of the wave profile £ is the remarkably simple expression 


_ 1 + n 2 B 2 
n 2 B y 


(40) 


The type-(l,l) tensor B cr u = B \£=£ crit has a number of interesting properties. Firstly, dC, is an eigenform 
of B 2 rlt , and it follows that (1341) . (1551) . (1551) are the only independent algebraic conditions on C that arise from 
the above analysis. Secondly, the bound (1551) is saturated by B = J-'crit) and the significance of this result is 
revealed by recalling that the term within the square root in (EH) may be written as a determinant : 

1 - k 2 X - k a Y 2 /4: = det (/ + kB) (41) 


where I is the type-(l, 1) identity tensor. Thus, 


det (I + nB crit) — 0 (42) 

and, hence, — 1/re is an eigenvalue of .Fcrit- Morever, det(I + kB) = det (I — kB) because B is generated from a 
totally antisymmetric tensor (the 2-form F ), and it follows that 1 /k is also an eigenvalue of J r cr it. Furthermore, 
the eigenvectors W + , W- satisfying 

B cri ,W+ = -W+, 

B clit W. = --W-, 

K 

are null with respect to the spacetime metric, 

W+W+ = 0, WlW- = 0 (45) 


(43) 

(44) 


which follows because nW + B cv itW + = W + W + , KW-B cvl tW- = — W-W- and UBU = iijiuF = 0 for all U. 
Finally, using (1501) . (1511) . (1501) it can be shown that 

dCBi it = XdC (46) 


where 

BKX + ^B 2 ) 

k 2 B' 2 (1 + k 2 B 2 ) ' 

However, (1551) . (1551) yield dC,B 2 rlt W + = dCyV+jn 2 , dC,BB t W- = d(W-/n 2 , respectively, and therefore 


(47) 


d(W+ = 0, dC,W- = 0 (48) 

since 1/n 2 . 

The above considerations demonstrate that the pair {W+, W-} is a natural basis for constructing solutions 
to (1551) . (1551) . (1551) . (1551) . We can choose W + W- = —1/2 and choose W + ,W- to be future directed without loss 
of generality and, using the above results, it is easy to see that 


C = exp 


f — r) W+ + exp ( 
\rriK J \ 


—rW_ 

tyik ) 


(49) 


satisfies (1341) . (1351) . (1361) . Furthermore, E3D is satisfied because VpW+ = 0, V^,kF_ = 0. The non-perturbative 
nature of (1491) in k is clearly visible in the argument of the exponentials. 

The asymptotic behaviour, as r —> ±oo, of a test electron with 4-velocity EH is determined by the pair 
{W +: W-}. The 3-velocity u = (u x , u y ,u z ) of the electron is related to the 4-velocity C as 


C = 


1 (d d d d \ 

71 =^ [at + Ux di + Uy d~ y + Uz Yz) 


(50) 
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where u 2 = u 2 + u 2 + u 2 z . Hence, recalling q < 0, it follows that 

dxW- 


lim u = 


v 

dxW- 


+ 


dtW. 


lim u = 


+ 


dtW- 


( 51 ) 


where dxW + = {dxW + . dyW + , dzW + ), dxW- = ( dxW _, dyW-, dzW-). Inspection of (1151) . C5l) shows that, 
without loss of generality, the 4-vectors W + /dtW+ and W_ /dtW- may be parameterised as 


W+ 

dtW+ 


u u r - - 

~pr, + v ^~ + Vl -« 2 
at oz 


cos a_t_ ——|- sin a+ — 
ox oy 


W- 

dtW- 


d d 
dt +V d^ 


VI 


- d d \ 

— v z cos a_ ——h sm a- — 

\ ox oyJ 


(52) 

(53) 


where a+,a_ are the angles (from the s-axis) of the incoming and outgoing trajectories, respectively, projected 
into the x — y plane. Consideration of (l43l) . dill) . (l52l) . (l53l) yields 


sin(a± — ip) = ± 


v/1 + k 2 B 2 cos 2 9 


cos(a± — (p) = 


kB cos 9 


Vl + k 2 B 2 cos 2 9 


(54) 


where 9 , ip are spherical polar angles that specify the orientation of the ambient magnetic field relative to the 
z- axis: 


B x = B sin 9 cos <p, B y = B s\n9 suup, B z = Bcos9. (55) 

The above results may be used to qualitatively estimate the behaviour of an electron captured in the peak of a 
pulse propagating at speed v. We see that the electron emanates from the core of the pulse with 3-velocity = 
(Vi — v 2 cosa_, VI - v 2 sin a- , v ) where a_ = p + arccot(— kB cos 9). Hence, for physically reasonable values 
of nB 1 the electron is strongly accelerated to near the speed c = 1 and ejected from the pulse at a finite angle 
to the direction of propagation of the pulse. In particular, «_ = (— sin(9 sin c/j, kB sin9cosip, 1) + 0{n 2 B 2 ) 
and the electron is ejected from the core of the pulse at the azimuthal angle p + 7t/2 + 0{kB). 


4 Conclusion 

Two distinct non-perturbative results have been presented that address particle acceleration in non-linear elec¬ 
trodynamics. In the absence of an established theory of radiation reaction in the context of non-linear electrody¬ 
namics, we focussed our attention on a simple matter model compatible with stress-energy-momentum balance 
and explored test particle motion in that context. In Section [2] we obtained an expression for the maximum 
energy gained by a test electron in a non-linear density wave in a magnetised plasma. The expression is valid 
for a wide range of physically permissible theories of non-linear electrodynamics encoded by a Lagrangian of 
the two electromagnetic invariants (including Born-Infeld electrodynamics). However, if Born-Infeld electrody¬ 
namics emerges at sufficiently strong field intensities then, as shown in Section [3j an outcome may occur that 
cannot be realised using perturbative considerations based the linear vacuum Maxwell equations. Under the 
appropriate conditions, a test electron can ‘surf’ a critically intense Born-Infeld electromagnetic plane wave and 
be strongly accelerated by the wave. 
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